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Abstract 

We derive a upper bound on the free energy of a Bose gas system 
at density g and temperature T. In combination with the lower bound 
derived previously by Seiringer |15j . our result proves that in the low 
density limit, i.e., when a^g 1, where a denotes the scattering length 
of the pair-interaction potential, the leading term of A/ the free energy 
difference per volume between interacting and ideal Bose gases is equal 
to 47ra(2g^ — [g — Qc]+)- Here, gc{T) denotes the critical density for 
Bose-Einstein condensation (for the ideal gas), and [•]+ = max{-,0} 
denotes the positive part. 



1 Introduction 



The ground state energy, free energy are the fundamental properties of a 
quantum system and they have been intensively studied since the inven- 
tion of the quantum mechanics. The recent progresses in experiments for 
the Bose-Einstein condensation, especially the achievement of Bose-Einstein 
condensation in dilute gases of alkali atoms in 1995 [1], have inspired re- 
examination of the theoretic foundation concerning the Bose system, e.g., 
[H, [m, [13], [Z], m, [13, M and [Mj on ground state energy and [IS] on 
free energy. 

In the low density limit, the leading term of the ground state energy 
per volume was identified rigorously by Dyson (upper bound) [i3] and Lieb- 
Yngvason (lower bound) [13] to be Anag^, where a is the scattering length 
of the two-body potential and g is the density. We note that Anag'^ is 
also the first leading term of AE the ground state energy difference per 
volume between interacting and ideal Bose gases. (The ground state energy 
per volume of ideal Bose gas is zero). 

On the other hand, the first leading term of A/ the free energy difference 
between interacting and ideal Bose gases is the second leading order of the 
free energy. More specifically, when a^^ <C 1, where a denotes the scattering 
length of the pair-interaction potential, then 

f{Q, T) = fo{g, T) + 4^(2^,2 -[g- g^]l) + oiag") (1.1) 

Here, / is the free energy per volume of interacting Bose gas, /o is the one 
of ideal Bose gas and gc{T) denotes the critical density for Bose-Einstein 
condensation (for the ideal gas), and [•]+ = max{-,0} denotes the positive 
part. The lower bound on / has been proved in Seiringer's work [15] . In 
this paper, we prove the upper bound on / and obtain the main result (jl.ip 
The trial state we use in this proof is a new type, which was first used 
in [M]. Let (po be the ground state of ideal Bose gas system. The trial state 
(pure state) Yau and Yin constructed for interacting Bose gases in [16] is 
almost equal to the following one 

k vr^^/g 

+ ^ Cfc4al;^ao Cg |(/>o), (1-2) 

k 

(with suitably chosen c and A). This trial state (pure state) in [16] is used 
to rigorously prove the upper bound of the second order correction to the 
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ground state energy, which was first computed by Lee- Yang [9] (see also 
Lee-Huang- Yang [8] and the recent paper by Yang [T7] for results in other 
dimensions. Another derivation was later given by Lieb p!0] using a self- 
consistent closure assumption for the hierarchy of correlation functions.) 
We can rewrite the pure state (|1.2p as follows 



(II21) = %o)%v^)l0o) 



where 



(0,0) 



(0,v^) 



exp 



exp 



E 



Cfc4a;Lfeao oo 



(1-3) 



(L4) 



Y Yj '^\J ^k+v/2^-k+v/2a\+v/2^'^-k+v/2'^v ^0 



We note: -F'(o,o) represents the interactions between condensate and conden- 
sate, since two particles with momenta zero are annihilated (og Cq ) and two 
particles with high momentum are created (a|.ol^). Similarly P(o,^) repre- 
sents the interaction between condensate and the particles with momentum 
~ 0{q^^'^), since in this operator one particle with momentum zero and one 
with momentum ~ 0{q^^'^) are annihilated (a^ ag ) and other two particle 
with high momenta are created. 

In this paper, we construct a trial state of the similar form. More specif- 
ically, let L/ be Gibbs state of ideal Bose gas at temperature T, the trial 
state we are going to use is very close to the following one 



r ~ ( -P(gl/3^gl/3)-P(0^gl/3)-P(0,0) j T/ ( ^'(gl/3,gl/3)-P(g^gl/3)-P(0,0) 



(L5) 



where 



Pi 



(0,0) 



-P(0,el/3) 



-P(£)l/3,gl/3) 



exp 



exp 



exp 



Y Cfc4o-feOo ao 



L k 



(1-6) 



Y Y '^J^k+v/2>^-k+v/2a\+v 



t J 



/2^~k+v/2^v "-O 



Y Y ^ ^k+■!^^-k+^-k+^■^-k+^ 
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where constant 2 comes from the ordering of . As one can see -P(o,o) rep- 
resents the interactions between condensate and condensate, P^q ^1/3^ repre- 
sents the interaction between condensate and tlic particles with momentum 
~ 0{q^/^) and -P(^i/3 ^1/3) represents the interaction between the particles 
with momentum ^ 0{q^^^). 



2 Model and Main results 

2.1 Hamiltonian and Notations 

Wc consider a Bose gas system which is composed of N same bosons and 
confined in a cubic box A of side L. The Hilbert space 'Hn,a for the system 
is the set of symmetric functions in L^(A^). The Hamiltonian is given as 

N 

HM,A = -J2^i+ E (2-1) 

i=l l<i^j<N 

Here the two body interaction is given by a smooth, symmetric non-negative 
function V{x) of fast decay. In particular, it has a finite scattering length, 
which wc denote by a. As usually, we denote by ^ a) Hamil- 
tonians with periodic (Dirichlet) boundary conditions. 



'27^r7^3 

1 apeii^c ui iv io iv . — \ 

^ we have 

J_ T F(p) = — T F(p) / 

The Fourier transform is defined as 



The dual space of A is A* := (^Z)'^. For a continuous function F on 



and then 



1^1 peA 



where Jrs (x) is the usual continuum delta function and the function d\* (p) = 
|A| = L*^ if p = (otherwise it is zero) is the lattice delta-function. We will 
neglect the subscript, the argument indicates whether it is the momentum 
or position space delta function. In general we will also neglect the hat in 
the Fourier transform. To avoid confusion, we follow the convention that 
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the variables x, y, z etc denote position space, the variables p, q, k, u, v etc. 
denote momentum space. We also simplify the notation 

i.e. momentum summation is always on the A*, it will be more convenient 
to redefine the bosonic operators as 

1 t 1 t 

(without changing the notation) i.e. from now on we assume that 

r ti- t_t _/l Sip = q 
[ap , flgj - % GgUp - I Q otherwise. 

Thus our Hamiltonian in the Fork space = ®n'Hn,a, is given by 

Ha = ^p'^alup + XI (2-2) 

2.2 Free energy 

The free energy per unit volume of the system at temperature T = f}''^ > 0, 
density q = A^/|A| > in cubic box A is given by 

/(^,A,/3) = -^\n{Txn^^^Exp{-PHN,K)) , (2.3) 

Let f^{g,A,l3) and f^{Q,A,l3) denote the free energy per unit volume of 
the system with periodic or Dirichlet boundary conditions. Furthermore, 
we denote by f{Q,P) the free energy (per unit volume) in thermodynamic 
limit, i.e., |A|, iV — > oo with g = iV/|A| fixed, i.e., 

f(^\g,P)^ lim r(^)(^,A,/3) (2.4) 
|A|-^-oo 

As mentioned in the introduction, in this paper, we give a upper bound on 
the leading order correction of f{g,(3), compared with a ideal gas, in the 
case that a^g is small and lim^^o/?^^^^ € (0, oo). We note that a^g and 
Pg"^/^ are dimensionless quantities. 
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2.3 Ideal Bose gas in the Thermodynamic Limit 

In this section, we review some well known results. In the case of vanish- 
ing interaction potential (V = 0), the free energy per unit volume in the 
thermodynamic limit can be evaluated explicitly. Let C denote the Riemann 
zeta function. It is well known that when q^/'^(3 > (47r)~"^^(3/2)^/^, i.e., g is 
greater than critical density Qc, 

Q>Qc= (4^/3)-'/'C(3/2) (2.5) 

the free energy in the thermodynamic limit is given as 

C\e.P) = / ln(l - e-P^')d?p (2.6) 

On the other hand, when Q < Qc, 

fo'^''\o,P) = Qf^ + 7:^ [ HI - e-^'~P'-^^)d'p (2.7) 

Here fi{Q, /5) < is determined by 

Note: when g > Qc, fJ-{g, (3) is defined as zero. 

It is easy to see the scaling relation: fQ^^\Q,f3) = Q^^"^ fo^^\Q'^^^ 
and the ration Qc/q only depends on dimensionless quantity g'^^^P, i.e., 

gjg = (4^)-3/2<^(3/2)(^2/3^)-3/2 ^2.9) 
Let P{g) be a function of g, we define R[f3] as the ratio gc/Q in tlie limit 

R[(3] ^ lim g,if3)/g = lim(4^)-3/2^(3/2) (g^/^ p^g))"'^^ (2.10) 

2.4 Scattering length 

In this paper, we use the standard definition of scattering length, as in [H] . 
[2]) [I!, [E]) [S], [E], [15]. Let 1 — t;; be the zero energy scattering solution, 
i.e., 

- A{l-w) + V{l-w) = (2.11) 
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with < w < 1 and w{x) — > as |x| — > oo. Then the scattering length is 
given by the formula 



a :- 



1 

4tt 



V{x){l — w{x))dx 



With (I21T]) . we have, for p 0, 



[V{1-w)Up\ 



This implies that if V is smooth, then 



dwr, 



dp 



< const. (\p\ ^ + \p\ ^) 



(2.12) 



(2.13) 



(2.14) 



On the other hand, because V{1 — w) > 0, so for Vp, 

[Vil-w)l\ < I Vil-w). 

Then with (|2.12p . i.e., J V{1 — w) is equal to 47ra, we obtain the following 
bound on Wn 



\wp\ < 47ra|p|' 



(2.15) 



2.5 Main results 



THEOREM 1. In the temperature region where linig^Q g'^^^ P{q) S (0,cxd) 
and in thermodynamic limit, for fixed scattering length a, we have the follow- 
ing upper bound on the free energy difference per volume between interacting 
Bose gas f^ig^fi) and ideal Bose gas fQ'{g,(3), 



lim,_o {/""{q, P) - fois, /?)) 0'^ < 4™(2 - [1 - R[p] (2.16) 

where R[P] is defined in (I2.10p as the ratio gd Q in the limit g ^ 0. 

It is well known that the effect of boundary condition in the thermody- 
namic limit is negligible, i.e.. 



Mq,P) ^ foiQ^P) = fo{e,P) and f{g,(3) ee f^{g,P) = f''{g,(]) (2.17) 

So with the lower bound on f^{g,(3) in Seiringer's work in [15], we can 
obtain the following result. 

COROLLARY 1. With the assumption of TheoremUl we have: 

lim {f^{g, (3) - fl^ig, (3)) g-^ = 4™(2 - [1 - ]^), (2.18) 
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3 Basic strategy 

3.1 Reduction to Small Torus with Periodic Boundary Con- 
ditions 

To obtain the upper bound to the free energy, we can use the variational 
principle, which states that, for any state r^(^) {Hn Hn) satisfying 
Dirichlet (Periodic) boundary condition, the following inequality holds. 

/^^^^(^, A, /3) < ^Tr^.,A ^^,Ar^(^) - \^pS{T^^''^) (3.1) 

Here, SiV) = — TrFlnF denotes the von-Neumann entropy. Hence, to prove 
Theorem [H one only needs to construct a trial states T^{q, A, /?) satisfying 
Dirichlet boundary condition and the following inequality: 

< 47ra(2 - [1 - ]^) (3.2) 

Furthermore, the proper trial states in thermodynamic limit (A oo) 
can be constructed by duplicating the proper trial states in the small box 
(A = Q~'^, c > 0) with Dirichlet boundary condition. Hence, the following 
Proposition one implies our main result Theorem one. 

Proposition 1. In the temperature region where liuig^Q g'^^^ P(q) E (0,oo), 
for fixed scattering length a, there exist A, |A| > and trial states 

T^{q, a, (3) satisfying Dirichlet boundary condition and the following in- 
equality, (Here N = \A\g) 

^,^0 (^Tr/7jv,Ar^ - jl^S{T^) - /?(f^,/3)) g'^ 
< 4Tra{2-[l-R{(3)]l), (3.3) 
where R{I3) is defined in ()2.10p . 

On the other hand, the next lemma shows that a Dirichlet boundary 
condition trial state with correct free energy can be obtained from a periodic 
one. 

Lemma 1. Let the volume |A| be equal to g"^^^"^^ . In temperature region of 
theorem one, if 

f^'ig, A, P) < const, g^^^^ (3.4) 
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then for revised box A* and density g* : 

|A*| ^ |A|(1 + 2/1/^20)3^ ^* ^ ^ 2q^^I^^^)-\ (3.5) 
we have f^{Q*, A*, /3) bounded from above as follows 

IiE,^o {f'^ig*, A*, /?) - A,/?)) ^,-2 < (3.6) 

We note: |A*| > {q*)~^^/^'^^ . The construction of a periodic trial state 
yielding the correct free energy upper bound is the core of this paper. 
We state it as the following theorem, which gives the upper bound on 
f^ig, A, P) in ^ and M . 

THEOREM 2. Assume lim^^o e-^^^^ G (0,oo) . For |A| = g'^^/'^^, N = 
\A\g, there exists trial state T{g,A,f3) satisfying 

m,^o (^Tr H^T - j^S{T) - f^ig, /?)) g-' < 4na{2 - [1 - R[(3] ]l) 

(3.7) 

It implies 

m,^o WiQ, A, /3) - f^{g, /5)) < 4^„(2 _ [i _ ]^) (3.3) 
3.2 Proof of Proposition 1 

With Lemma [T] and Theorem [2l we can prove Proposition 1 as follows. 

Proof. Using the temperature function (3 in the assumption of Proposition 
1, we define a new temperature function (3 as follows 

Aq)=I3{q*), (3.9) 
where g* = g{l + 2^4i/i20)-3^ ^ -^^ 

Insert the result in Theorem [2] into Lemma [TJ With the definition of A* , 
g* in Lemma [T113.5p . we obtain in the inverse temperature f3{g), 

iii^^^o [f'^iQ*, A*, ^) - fP{g, P)) g~^ < 47ra(2 - [1 - R[^] ]\). (3.10) 

Using the well known results of ideal Bose gases: 

fP{g, P) = f^{g, P) = f^{g\ M + oig"/^)), (3.11) 

we can replace /o'(f?, P) in (|3.10p with fQ{g*, P), i.e., 

^e^o [f'^iQ*, A*,P)- foi6*, P)) < 4^a(2 - [1 - R[P] (3.12) 
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Then by we obtain R[l3] = R[f3], so 

TSI,^o {f''{Q*, A*, Piel) - f^{Q\ P{q*))) < 47ra(2 - [1 - R[P] ]^) 

= 4Tra{2 - [1 - R[p]]l) 

41 

At last, with the facts: A* > (^*)~2o and the hmit ^ — > is equivalent to 
the limit g* — > 0, we arrive at the desired result (|3.3|) . ■ 

3.3 Reduction to Pure States 

The Lemma [T] can be proved with standard method as in [16j and we leave 
the proof in section [12.11 In this subsection, we introduce the basic strategy 
of proving Theorem [21 With the assumption of Theorem [21 we have 

A=[0,Lf,L = g-^,N = Q-^and \im g'^/'^p £ {0,oo). (3.13) 

We first identify four regions in the momentum space A* which are rele- 
vant to the construction of the trial state: Pq for the condensate, Pl for the 
low momenta, which are of the order g^^^; Ph for momenta of order one, 
and Pi the region between Pq and Pl- 

DEFINITION 1. Definitions of Po, Pi, Pl and Ph 

Define four subsets of momentum space A* = (27rL~^Z)^; Pq, Pi, Pl 
and Ph as follows. 

Po ={p = 0} 

Pi = {p G A*|0 < IpI < eL^^/^} 

Pl ={pe A*\eLg^^^ < \p\ < vI'q^^^} (3.14) 

Ph ^{peA*\eH<\p\<VH} ^ 
where the parameters are chosen as follows 

£L,'nL,£H,'nH = and r] = 1/200 (3.15) 

We remark that the momenta between Pl and Ph are irrelevant to our 
construction. 

DEFINITION 2. Definition of M, M and 

Let P denote Pq U Pl U P/ U Ph- We define M as the set of all functions 
a : P ^ N U such that 

Y^a{k)=N (3.16) 

fceP 
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For any a G M, denote hy \a) G 'Hn,\ the unique state (in this case, an 
N -particle wave function) defined by the map a 

k&P 

where the positive constant C is chosen so that \a) is L2 normalized. 
Moreover, we define M as the following subset of M 

M = {aG M|supp (a) C PqUPiUPl and a{k) < rric for V/c G Pl}, 

(3.17) 

where rric is defined as 

mc ^ e-'"^ = (3.18) 

Clearly, we have 

alak\a) = a{k)\a), yk e P (3.19) 
With the definition of Pq, Pl, Pj and Ph, for a G M, we define G M as 

iVa = a(0)a(0) + ^ 2a{u)a{v) (3.20) 

To prove Theorem [21 first, we show that, with q's in M (I3.17p . we can 
construct a trial state Fq satisfying (j3.7p . but with wrong coefficient in RHS. 

Lemma 2. For A = [0, LJ'^, L = q eo ^ N = q 20 and lim^^o P £ 
(0,c«). There exists a state Tq(q, [3) having the form: 

To = 5a(f?,/3)|a)(a|, ffa(^,/3) G M (3.21) 
and satisfying (It is Vq = /y(x)(ia;^, noi 47ra m i/ie r/is.j 

Ik^^,^o (^Tr F^vTo - ^'^(ro) - /o(^,/3)) Q'^ < ^o(2 - [1 - ]\) 

(3.22) 

Furthermore, with defined in (j3.20p , ^/le coefficient function g^ satisfies 
lim iV-^A^aSa = 2 - [1 - i?(/3) ]^ (3.23) 

We remark: actually Tq is very close to Tj the canonical Gibbs state of 
ideal Bose gases. The state ro(^, /3) satisfies ()3.22p . but for most potential 
V{x), Vq = JV{x)dx'^ is strictly larger than 47ra. So we need to improve 
Tq. To do that, we need to replace |a)'s in M with some ^^'s. To introduce 
we start with dividing the momentum space as follows. 
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DEFINITION 3. Definitions of Bh{u), Bl{u) 

Let >cl,>ch > 0. Divide Pl and Ph (|3.14p into small boxes (could he 
non-rectangular box) s.t. the sides of the boxes are about and q^'^ . We 
denote the box containing u as Bh{u) when u S Ph or Bl(u) when u £ Pl- 

DEFINITION 4. Definition ofMa _ 

For any a G M, we define Ma as the set of the f3's in M (Def. that 

1. Ifke Pi, then p{k) = a{k). 

2. There is at most one k in each Bl or Bh satisfying (3{k) ^ a{k). 

3. If I3{k) / a{k), then 

p{k) = a{k) - 1, for k e Pl (3-24) 
(3{k) = a{k) + 1 = 1, ior kGPn 

For each a £ M, we will construct a normalized pure state which is 
linear combination of /? S M^, i.e., 

= Yl E = 1 (3-25) 

To prove Theorem [2l i.e., to improve the Tq in Lemma [2j we choose the 
correct trial state T as follows: 

r = Y ga\^a){^a\, (3.26) 

where we choose g^ in (j3.2ip and in (|3.25p . 

With proper x/, and x/^, AS the entropy difference between Tq in (j3.2ip 
and r in p.26p can be proved to be much less than Ag'^. 

Lemma 3. Let A = q~'^^/'^'^, kl < 5/9 and kh < 2/9. Then for any 
{^'a,a G M} having the form (j3.25p . we have 

m,^o [ - S{T) - (-S(ro))] {Ag^r' = (3.27) 

We remark: the assumptions < 5/9 and < 2/9 implies: 

In the next theorem, we will show, for each a G M, there exists a pure state 
^a- Comparing with |a), the new pure state l^'a) lowers the total energy 
by about (Vq — 4:TTa)NaA^^ . The construction of the pure state yielding the 
correct total energy is the core of the proof for theorem [2l 
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THEOREM 3. Let 1/2 > > 4/9 and kr > 1/9. For any ae M, there 
exists having the form ()3.25p and satisfying: 

{'^a\HN\'^a) - {a\HN\a) + (Vb - 47ra)iVaA""^ < e^^^A 
where the Eg is independent of a and limg_^Qeg = 0. 

3.4 Proof of Theorem [2] 

Proof Let 1/2 > > 4/9 and 2/9 > >ch > 1/9. We choose trial state T 
(|3.26|) with ga in Lemma [2] (|3.21|) and ^^'s in Theorem [31 Then Combine 
Theorem [3l Lemma [3] and Lemma [2l ■ 

This paper is organized as follows: In Section HI we rigorously define 
^'q's and the trial state T. In Section 5, we outline the Lemmas needed 
to prove Theorem [3l In Section 6, we estimate the number of particles 
in the condensate and various momentum regimes. These estimates are 
the building blocks for all other estimates later on. The kinetic energy is 
estimated in Section 7 and the potential energy is estimated in Section 8-11. 
Finally in Section 12, we prove Lemma [U [2l [3l 



4 Definition of the trial pure states ^^'s 

In this section, we give a formal definition of the trial pure state ^'^'s for 
Theorem O We start with defining an (pair creation) operator Ap'g: 

Ap;^ -.M^M, u,v e PqU Pl, p,q e Ph and u + v = p + q (4.1) 

For l3 eM, 

\A;fg(3) = Calalaua^lp), (4.2) 

where C is positive normalization constant. The operator Aplq annihilates 
two particle with momenta in Pl or Pq and creates two particles with mo- 
menta in Ph- We note: the total momentum is conserved. 

For simplicity, the pure trial state ^'^ will be of the form J^fs^Ma faiP)\P) 
where fa is supported in C Ma which we now define. 

DEFINITION 5. Definition of nontrivial subset in Pl and Ma 
Let A he subset of Pi, it is called non-trivial when 

1. If ui, U2 £ A, Ui 7^ Uj (i 7^ j), then ui + U2 ^ 
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2. If ui, U2, U3 £ A, Ui ^ Uj(i 7^ j), then ui + U2 ^ 

3. If ui, U2, U3, £ A, Ui ^ Uj (i 7^ j), then ui + U2 u^ + u^. 

Then recall Ma in Def. 0. For a € M , we define the subset C Ma 
as the smallest set with the following properties. 

1. Denote the following subset of Pl as ^^(7,0), 

Pl(7, a) = {n G Pl : 7(n) < a{u)}. (4.3) 
Pl{'y,a) is non-trivial, for any 7 G 

2. a£ Ma 

3. If 13 £ Ma and 7 = A^^^_pf5 G Ma, then 'y £ Ma- 

4. Ifpe Ma, 7 = G Ma and 

(a) PL{'~f,a) is non-trivial 

(b) Pi-p) = a{-p), Pi-q) = a{-q) 

then 7 G Ma- 

Note: The set Ma is unique since the intersection of two such sets Ma 1 
and Ma,2 satisfies all four conditions. The properties of non-trivial set have 
no physical meaning, but they can simplify our proof and calculation. 

We collect a few obvious properties of the elements in Ma into the next 
lemma. 

Lemma 4. By the definition of Ma, any (3 G Ma has the following form: 

m n 

P = X{^^Z^X{^,-P.-^ (4-4) 

i=i j=i 

where Ui G Pl^Pq, ki G Ph for i = 1, - ■ ■ , 2m and pj G Ph for j = !,••• , n. 
And 

Pi / ki / ±kj for i j and ki / ±pj for Vi, j (4-5) 

On the other hand, if{ui, (i = 1, • • • , 2m)}nPL is a non-trivial set of Pl, any 
[3 with form (|4.4|) and (|4.5|) belongs to Ma- Furthermore, one can change 
the order of the A's in (|4.4p . With the fact that the subset of non-trivial 



14 



subset of Pl is still non-trivial, we can see, if (3 belongs to Ma and has the 
form ()4.4p and ()4.5p . then 

Here A, B are any subsets of {1, • • • , m} and {1, • • • , n} 

DEFINITION 6. The Pure Trial State 

Recall function {1 — w) is the zero energy scattering solution of potential 
V , as in (|2.1ip . Define the pure trial state ols 

|*a) = E (4.7) 

/3eM„ 

where the coefficient 's are given by 



l3{k)>0 I3(k)>p{~k) f3{u)<a{u) 



^"(«=^"Vw n ^ n n '-'-'^ 

V kePH kePa uePL » i i 

Here we follow the convention y/x = \/\x\i for x < Q. For convenience, we 
define /(/?) = for (3 ^ M^. The constant Ca is chosen so that ^'q is L2 
normalized, i.e., 

{^^\^a) = 1, i.e., Yl = 1 

With choosing as above, we can obtain a few obvious identities of 
as foUows. 

Lemma 5. 1. If k £ Ph and f3 G Ma,^°'° G Ma, then 



2. If ui,U2 G Pl, U2 = ±tii or U2 G Bl{ui), ki,k2 G Pff and /? G Mq,, 
tften 7 = All'^lP i Ma, i.e., fa{l) = 0. 

3. Ifui,U2 G Pl U Po and U2 / ±ni, ki,k2 £ Pr , fi e Ma and Al\'^ll3 G 
Mq,, i/iera when fi{—p) = a{—p) and (3{—q) = a{—q), we have 



n^ZTJ) = %p/(/3) (4-10) 
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when l3{—p) / oi{—p) or l3{—q) / a{—q), we have 



< 



|A| V 1^1 



(4.11) 



Again the result 2 in Lemma [S] has no physical meaning, but it can 
simplify our proof. 

At last we give a brief explanation why we mentioned in introduction 
that r the trial state is very close to ()1.5p . With the definition of ^^^'^^ ' 
can see that almost equals to 



where 



Pi 



(0,0) 



-^(0,^1/3) 



exp 



exp 



exp 



*a ~ ^{gi/3,gl/3)-P(0,gi/3)^'{0,0)|a) 



L k 



(4.12) 



(4.13) 



E E - 



o„ a,. 



On the other hand the trial state T has the from ^aga\^a){^a\ and Fq 
5ci|a)(a|, so the trial state T is very close to the following one: 



r ~ ( -f'(^>l/^^?l''3)-^(o,^?l''3)-^(0:0) ) ^o ( -f(f)i/3,pi/3)-f(o,ei/3)-P(o,o) 



(4.14) 



As mentioned before, actually Tq is very close to Tj the Gibbs state of ideal 
Bose gases, so we obtain 



r ~ ( -Pol/3 ni/3-Poni/3-f'o,0 1 T/ ( ^„i/3 „l/3-Po,rtl/3-Po,0 



t 



(4.15) 



5 Proof of Theorem [3] 

Proof. Our goal is to prove 



{^a\HN\'^a) - {a\HN\a) + (Vb - 47ra)7V„ < e^^i^A 



(5.1) 
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First we decompose the Hamiltonian H^. By the rule 1 of the definition 
of M„, j3{k) is equal to a{k) for any k ^ Pj and /3 G Ma C Ma- So if 

ki € P/, /3, 7 G Mq, and (/3|afc^a^2'^fe'^^4l^) ^' ^^^^ '^^^ of /C3 and k4 must 
be equal to ki. On the other hand, since the particles with momenta in Ph 
are created by pair, the total number of the particles with momenta in Ph 
is always even. With these two results and momentum conservation , we 
can decompose the expectation value {'^a\HN\'^a) as follows: 

N 
i=l 

(5.2) 

where 

1. Habab is the part of interaction that annihilates two particles and cre- 
ates the same two particles, i.e., 

(5.3) 

2. Hj^j^ is the interaction between particles with momenta in P^: 

Pi = PoUPl (5.4) 

and 

where ui 7^ U3 or 144. 

3. Hj^jj is the part of interaction that involves two particles with mo- 
menta in Pj^ and two particles with momenta in Ph i.e., 

^LH= l^r^ X K<i-fei4i42«feiafe2 +C-C (5.6) 

ui,u2ePj^,ki,k2ePH 

+ X ^{yui-u2 + Vui-k2)ai^al^au2ak2, 

where ui ^ U2- 

4. Hhh is the part of interaction between particles with momenta in Ph, 

Hhh = |Ar^ X ^fc3-fei4i42^fe3«fe4> (5-7) 

where ki 7^ or ^4. 
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With these definitions, we can rewrite the total energy of \a) as: 



N 



{a\HN\a) = {a\ ^ -Aj|a) + {a\Habab\a) (5.8) 



i=l 



The estimates for the energies of these components of in ()5.2p about ai^e 
stated as the following lemmas, which will be proved in later sections. 

Lemma 6. The total kinetic energy is bounded from above by 

/^-aA -/^-aA -\\Vw\\lN^\A\-' <eig^A, (5.9) 

\i=l / \i=l fa 

where ei is independent of g and limg^o^i = 0- 

Lemma 7. The expectation value of H^bab is bounded above by, 

Lemma 8. The expectation value of Hjj^ is bounded above by, 

{H-Li)^^ < e'"'^ (5.11) 

Lemma 9. The expectation value of Hj^^ is bounded above by, 

{^lh)^^ + N^\^\~^pVw\\i < e2Q^A, (5.12) 

where 82 is independent of a and limg^o^2 = 0. 

Lemma 10. The expectation value of Hhh is bounded above by, 

(Hhh)^^ - N^\A\-^\\Vw% < esQ^A, (5.13) 

where £3 is independent of a and limg_^o^3 = 0. 

On the other hand, by definition of w in (j2.1ip and (j2.12p , we have 

\\Vw\\l - \\Vw\\i + WVw'^Wi = 0, Vo- \\Vw\\i = iira (5.14) 

Together with ()5.8p and ()5.9p - (j5.13p . we arrive at the desired result (jS.ip . 
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6 Estimates on the Numbers of Particles 

The first step to prove the Lemma [6] to Lemma [10] is to estimate the particle 
number of in the condensate, Pl,Pi, and Ph- This is the main task of 
this section and we start with the following notations. 

DEFINITION 7. Suppose Ui £ P for i = 1, . . . s. The expectation of the 
product of particle numbers with momenta ui, ■ ■ ■ Ug'. 

Q,(ni,n2,--- ,ns) = /n<a«.) = 11/3(^0 1 I ' (6-1) 

\i=l I f3GMa i=l 

DEFINITION 8. We denote by Ma{u) C M„ the set of (5's satisfying 
fj{u) = a{u), i.e. 

Mo,{u) = {P £ Ma : Piu) = a{u)} (6.2) 

Furthermore, with the definition of and Bh, we define M^{u) C Ma{u) 
as the intersection of Ma{v) 's of all v G Bl{u) or Bh{u), i.e., 

M^{u) = n^eBL{u)orBH(u)Ma{v) (6.3) 

We can see 

(3 £ M^{u) ^ P{v) = a{v) for Vt; G Bl{u) or Bniu) (6.4) 

Using (j3.24p . we can see that, for u £ Pl, the expectation of the particle 
number Qa{u) = {aia^)^^ is equal to a(^^) - E/3^a/<,(«) For k £ Ph, 

the expectation value of particle number Qa{k) is equal to ^^f^^Maiu) l/(/3)P- 
The following theorem provides the main estimates on Qa{u) and Qa{k). 

Lemma 11. For small enough g, Qa{u) and Qa{k) can be estimated as 
follows (u,ui,U2 £ Pl o,nd k £ Ph) 

< a{u) - Qaiu) = Yl l/(/3)l^ ^ co^st. g^-^"^, for u £ Pl(6.6) 

Furthermore, the probabilities of the combined cases are bounded as follows: 
(u, ui,U2 £ Pl and k £ Ph ) 

Y < const, g^-^"^ w\ienui^U2 (6.7) 

l3^Mc{ui)VjMc(u2) 

Y \fm' < const. g^-'^\wk\ (6.8) 
l3^Ma{u)UMa{k) 
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Proof. Proof of Lemma [TT] 

First, we prove (j6.5p concerning k G Ph- Using Lemma Hr ()4.4p - ()4.6p ). 
we have when (3{k) > 0, there exist 7 G Ma and u,v £ PlU Pq, p G Ph s.t. 

^^p7 = P and p = u + v — k (6-9) 

With the properties of /q, in Lemma [5lf (|4.9p - (|4. lip ). /«(/?) can be bounded 
as 

< 47H7(^;)A-2 \wkWp\ \fa{l)\'. (6.10) 

Then sum up /? ^ Ma{k), i.e., /3(A;) > 0, by summing up n, v and 7, we 
obtain: 

p^Ma{k) u,vePLUPo ■yeMa 

< 4:Q'^\wk\ maxjl^pl} (6-11) 

Then, with the result in (j2.15p : {wpl < 47ra|p|^^, we obtain 

Qa{k) < const. £>^~^''|wfc|, A; G -P/f (6.12) 

Using ()2.15p again, we obtain (j6.5p . 

Then, we prove (j6.6p concerning u G P^. Similarly, with Lemma 21 for 
any /3 ^ Mq,(u), i.e., P{u) = a{u) — 1, there exist 7 G Ma and v G Pl U Poj 
p,ke Ph s.t. dlJ]) holds. This implies (KIM . Using (f2T5]) and <^\k\, 
we have 

l^i'p'U^fcl < const. I A:|~^, when p. A; G Ph and |p + A;| ^ |A;| (6.13) 

Inserting this inequality and the bounds a{u) < rric = Q^^"^ into (I6.10p . we 
obtain: 

\fM < const. e-'^^lklSivWlfail)]' (6.14) 
Again, summing up /3, with 'y{v) < N, we obtain (j6.6p as follows 

^ |/a(/3)|'< ^ const. ^-3''|A;|-S(t-)A-Va(7)l'<£>'-'^ 

(6.15) 

Then, we prove ()6.7p concerning ui,U2 G Pl. For any /3 ^ Mq(ui) U 
Ma{u2), using Lemma m we can see that there are only two cases: 

1. there exist one 7 G Mq, pi,p2 G Ph and ^pi,'p2^7 = (3 
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2. there exist one 7 ^ Mq,(u2), f G Pl U Pq; ''^ / ^^2, Pi)i52 G -Pff and 

As before, with the properties of fa in Lemma [5] , the bounds on a(ti)'s 
{u £ Pl) and (j6.13|) . we have 

\Um' < const. Yl 0-'''\M-'\Ul)\' (6.16) 
+ const. Yl ^>-NW|ArVa(7)l' 

t^ePLUPo,7^Afc(M2) 

Using X]t)7(^) — ^ ^^"i (|6.6p . we obtain (j6.7p . 

At last, we prove (j6.8p concerning u £ Pl and k £ P^. For any /3 ^ 
Ma{u) U Ma{k), Using Lemma 01 we can see that there are only two cases: 

1. there exist 7 G M„, v £ PlU Pq, P & Ph and ^^'^7 = /3 

2. there exist 7 ^ Mo(n), ui, ^2 G -Pi U Pq, P & Ph and ^p|fc''''7 = /3 
Summing up v, p ot vi, V2, p, we obtain 

Y < const. Y y^7(^)7(^)A~^| (7)r 

l3^Mc{k)VjMc{u) v^PlUPo 7 

+ V 4Q^\wk\m3.x{\wp\}\Ul)\^ (6.17) 

— ' P&Ph 

With the result in (|2T3]) : |t(;p| < 47ra|p|-^ and E,;7(^) < 

we have: 

Y \fam' < const. j{u)0'-'^A-'\wk\+ Y V"''I^A.||/a(7)l' 

(6.18) 

At last using (j6.6p and the fact 7(ti) < a(ti) < g and A = 4i/20^ 
obtain the desired result (j6.8p ■ 

Moreover Qa{k){k £ Ph), has a more precise upper bound as follows. 

Lemma 12. For A; £ Ph, and Qaik) is bounded above by: 

Proof. First using Lemma [H we have that, for any /? ^ Ma{k), there are 
two cases: 
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1. there exist 7 G Mq, s.t., fc7 = /3 

2. there exist 7 G M^,, u ±v e Pl U Pq, p e Ph, s.t., ^p|^'''7 = p. 

So using the property of fa in Lemma [5l = Z]^^M„{fc) is 

bounded above by 

<q(0)2A-2«;2^ ^ 2a{u)a{v)A-^\wkWp\ 

(6.20) 

where p = u + v — k. Since ttip = w.p and |p + A;| < 2{q^/^^'^), with (j2.14p . 
we have 

\\wk\- \wp\ I < const, g^/^-'^'^ (6.21) 
Inserting this into (j6.20p . we obtain 

Qa{k) < NaA~^wl + //^-^"It^fcl (6.22) 

Combine with the fact l^^l < const. q~'^^, we obtain the desired result (j6.19p . 

■ 

At last, We collect a few obvious inequalities of fa into the following 
lemma. 

Lemma 13. .■ Recall the definition of Ma{k) or Ma{u) in Def. (16. 3p . 
the results (|6.6p and (|6.5p implies: 

Yl I' < i?'""'A£>'"" < q'^' for k G (6.23) 

and 

Y \fM\' < q'-^'^Aq^-^ = e"^ for uePL (6.24) 

Recall the definition of non-trivial subset of Pl in Def. 0. For any fixed 
u,v G Pl, define Ma''" C Main) n M„ (v) as follows 

M^^"" = {P G Maiu)nMaiv)\{u,vU PLiP,a)} is non - trivial} (6.25) 

Here PL{P,a) is defined in (|4.3p as the set {u £ Pl : f3{u) < a{u)}. Then 
using (j6.6p and ()6.7p . we /ia?;e 

E |/a(/3)|' < (6.26) 
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Similarly, for any fixed u ^ Pl, define M^'" C Ma{u) as follows 

M°''' = {(3e M„(u)|{n U is non - trivial} (6.27) 

Then using ()6.6p and (j6.7p . we have 



1/2 



(6.28) 



/asi, with (|6.5p anc? i/ie /aci a(0) — /?(0) < X^feeP^ (^i^)' /lawe (5a(0) 
and Qa{^, 0) hounded as follows 



and 



a(0) > Qa(0) > a(0) - ^^/^7V 
[a(0)]'>Q„(0,0)> [a(0)]'-iV2^5/6 



(6.29) 
(6.30) 



7 Proof of Lemma [6] 

In this section, we estimate the kinetic energy of by proving Lemma [6l 
Proof. By the definition, 

= E u^Qaiu) and = E ^^"(^) 

\i=l / u^PlUPiUPh \i=l /a uSPlUP/ 

(7.1) 

On the other hand Qa{u) < a{u), for u G Pj U Pl- So we obtain the LHS 
of (|5.9p bounded above by 

/^-aA -(y.-^) -\\VwfMA\-' (7.2) 

\i=l / <I-^ \i=l /a 

< Yl - l|Vu;||^iV„|A|-^ 

With the upper bound on Qa{k) in (16.19p . we have 



21) < iValAI-l 



Then using the fact hmg_+o 
plete the proof of Lemma [H 



\Vw\\l- E \M''k'\wk\' 

kaPn 

|2 V- IAI~U2| 



+ Q^^'^K (7.3) 



\^w\\i-Y.k^pJk\~^k'\wk\ 



0, we com- 
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8 Proof of Lemma [7] 

Proof. First we rewrite the expectation value of Habab 

{Habab)^a (8-1) 
l3eMc\ u uj^v J 

= |A|-i lvo{N'-N) + Y,Vu~vP{umv)] |/„(/3)|2 

On the other hand, 

{Habab) a = |A|-i I Vo{N^ ' + Vu-va[u)a[v) j (8.2) 

By the assumptions, Vy is positive for any \v\ <^ 1. Together with the fact 
P{u) < a{u) for n G Po U P/ U Pl, we have for P e Ma 

Vu-vP{u)P{v) < Vu-vaiu)a{v), when u, t; G Pq U P/ U Pl (8.3) 

Using this inequahty and the fact a{k) = for k £ Ph, we have 

{Habab)'!' a ~ {Habab) a 
- I^I^M Yl 2Vu-vQa{u,v) + Vu-vQa{u,v) 

\u^PH,vePH u,vePH 
Using the fact \Vu\ is no more than |Vo| for any u £ P, we obtain: 

{Habab)^^ - {Habab)a < 2VoQ Qa{v) < Q^'^^A (8.4) 



For the last inequality, we used (16. Sp . 

9 Proof of Lemma [8] 



We start the proof with the following identity for {au^au2(^u3a.u4)<S'. 
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Lemma 14. For any fixed momenta iii,2,3,4 o-n-d (3 € M^, define T{(3) to he 
the state 

\T{l3))^Cai^ai^au,auM3), (9.1) 

where C is the positive normalization constant when \T{f3)) ^ 0. Then we 
have 

(^U4 I ^U4^ \ P) 

(9.2) 

The map T depends on tii,2,3,4 and in principle it has to carry them as 
subscripts. We omit these subscripts since it will be clear from the context 
what they are. 

Proof. For any fixed tii,2,3,45 by the definition of ^'q, we have 



a\ai^ai^au^aui\^ a) = ^ /a(/?)/a(7)(7l4i42««3««4l/5) (9-3) 



By definition of Mq, one can see 

(7|ati42 0«3a«4l/3) 7^ ^ 7 = T{f5) 
Since \T{j3)) is normalized, the identity in Lemma O is obvious. 

9.1 Proof of Lemma [8] 

Proof. Using the fact \Vu\ < |Vo| for any u G R^, we can see 



niePj,ni7^U3,U4 



^Ml '^"3 '^"4 



We are going to prove: 

E 

UdPL 

E 

E 

Ui^Plj and ui^U3,U4 



'^Ul ^U2 '^"3 '^"4 



: 



(9.4) 



(9.5) 



(9.6) 
(9.7) 



(9.^ 
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First we note (|9.6p is trivial. Because if /5 G Mq,, then PL{(3,a) is non- 
trivial, which tells if f3{u) < a{u) then (3{—u) = a{—u). 

Then we prove (j9.7p concerning n2,3,4 G Pl- By definition of Mq,, one 
can see that if {I3\a\au^au:iaui\^) 7^ 0, then / M4 and 7 ^ Mq(m2)) i-e., 
7(112) < a(u2). Furthermore, with the definition of fa (|2.3p . we have 



'a(n3)a(u4) 



/ /3(0)a(u2) 
Combine with Lemma [HI we obtain 



/(7) 



Applying (j6.6p in Lemma W\\ we obtain 

ao'^u2^"3^"4/ < const. a(u3)a(ti4)£»^~^'', 



which implies ([9] 
At last, we prove 



Similarly, we have 



^u-i O^Ui 



7^MQ{ni)UM„(M2) 

Again, using Lemma [TTl we obtain 



< const. a{u^)a{ui)Q 



2-8r) 



which implies ()9.8p . At last, combine (j9.6p - (|9.8p and we obtain 

I {Hll). I < 9'"'^ 



(9.9) 



(9.10) 



(9.11) 



(9.12) 



(9.13) 



(9.14) 



10 Proof of Lemma [9] 

We start the proof with estimating {a^u^a^u^ak-i^ak^)-^^ when ui = ±U2 € Pl- 
By the definition of M^, if /? G Mq,, u & Pl and /3(n) < «(«), then /3(ti) = 
a{u) — 1 and f3{—u) = a{—u). So we have: 

(ai^ai^ak^ak^)^^ = 0, for V/ci, ^2 G -Pff, = ±U2 G -Pl (10.1) 
For the other cases, we leave the bounds in the following lemma. 
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Lemma 15. Recall Pj^ = PqL) Pl- For u,ui,U2 € P^ and k,ki,k2 G Ph, 
we have 



kePH 



[10.2) 



E 



U.l,2<^Pj,Ul^±U2 



U\,U2&Pr 



Uly^U2 

E 

Ul,U2&Pj^,ki,k2&PH 



(^Il4i"«2"fc2)*c 



:io.4) 



where £4, £5, o'^e independent of a and limg_^o = /or i = 4, 5, 6. 

Proof. Proof of Lemma [9] 

Combine the bounds in (fTIU]) . (fTn:2l) . (fTIOI) and (flHSD. ■ 

10.1 Proof of Lemma [T5] 

Proof. First we prove (jl0.2p concerning u £ P^ and A; G P^- By the defini- 
tion of Ma, if (atolM'^fe^-A:)*;! 7^ 0) then it must be zero. With the property 
of fa in Lemma [5l (j4.9p . one can see that if (/3|aQolQOfca-fc|7) / 0, then 



/a(7) = -WkV7{0?-l(.0)^-'UP) 
Together with Lemma [HI we have 



(10.5) 



(alalaka^k)^^ = -Wk Yl {P{Of - P{0)) A'^lfam^ (10.6) 



Recall the definitions of M^'s in Def. H One can see if /3(0) > 1, then 



P G Ma, e Ma is equivalent to ^ G M^(/c) n M^{-k). So, we have 

{alalaka.kh^ = -Wk (/3(0)2 - /?(0)) A" (10.7) 

/3eA/<?(A:)nAf^^(-fc) 
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Using the bound on Yl^gMB(k) \fa{P)\ (|6.23p and the bounds on QaiO), 
Qa(0,0) in (lOQ]) and We obtain that 

J2 (/?(0)2-/3(0))|/,(/3)|2-a(0)2 < 0(^,1/6^2) (10 8) 

/3eMS(fc)nMS(~fc) 

Insert (jlO.Sp into (jl0.7p . Then summing up /c € Ph, with u = 0, we obtain 



[10.9) 



-VkWkA-^ + \\Vw\\i 

kdPn 



0, we obtain 



Combine with the fact hm^^o 
the desired result (|lU.2p . 

Next, we prove (jlO.Sp concerning ui ^ ±n2 G Pj^ and ki,k2 G Ph- Using 
the result 2 in Lemma O one can see 

(aj^^aj^^afciflfca = when U2 G Bl{ui) (10.10) 
So from now on, we assume U2 ^ Bi{ui). Using the property of fa in 
LemmaO we have that if (/^lat^atjOfciafca It) 7^ and /?,7 G M^, then 

/(7) = c;3y^y^\/^M^M/(/3) (10.11) 

Here C/j depends on /3 and C/3 is no more than 2. Especially, when (3 G 
M^(— ^1) n Mq,(— ^2), C/3 = 2. Again with Lemma [TH for fixed ui, U2 ^ 
Bl{ui), ki and /c2j we have 

(alial2afciafc2)i'a = ^/^w^^/^w^ X] C'/3/3('"l)/?(^^2)|/(/5)|^ 

(10.12) 

First, using the facts \ki + ^2] < 2Q^^^r]L and the bound on dwp/dp ()2.14p . 
we obtain Itu^j — im^jI < then 

IV-W^fciV-W'fca +Wki \ < Q^^^ (10.13) 
Insert (jlO.lSp into (|10.12p . we have 

{alalak.aj,,)^^ = {-wj,, + ©(^^^Z^)) ^ C^/3(ni)/3(n2)|/„(/?)p. 

(10.14) 
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Now we bound 



^ C^(3iui)(3iu2)\fam^. 



/3eA^- /36A^- 



In the case P ^ Ma{—ki) H Ma{—k2), using the result in (|6.5p and |C^| < 2, 
we have 



J2 CpP{m)P{u2)\fam' 
/3^Af<,(fei)nMa(fc2) 



< ga{ui)a{u2) (10.15) 



In the case /3 e M^i-ki) D Mo,{-k2), we have = 2. Recall M^''""^ in 
([OS]) and (lOTl) . By the definition of M„(Def. [5] rule 4), we have 

Cpp{ui)p{u2)\um^ = 2p{ui)p{u2)\um\ 

(10.16) 

where A is the following subset of Ma- {pi,P2,P3,P4 ■= ui,U2,ki,k2) 



A = {r\UM^{p,)) n (nLiM«(-A:,)) nM«i' 



(10.17) 



Using the results in Lemma E] and Lemma [TSi; ([63]), (IHIH . (fHIMl) . (IHIM 
and (j6.28p ) and a(n) < rric for u £ Pl, we obtain that if ui,U2 G Pl 



^/?(ni)/3(n2)|/o 



and if ui = 0,U2 £ Pl, we have 



a{ui)a{u2) 



< 0(f)i/6-6'') 



(10.18) 



^ l3{ui)f3{u2)fa(.P) - a{ui)a{u2) 
l3eA 



(10.19) 



Inserting (110. 15p . (110. 18p and (110.191) into (110.141) . with the fact \wp\ < 
47ra|p|~^, we obtain that for ui,U2 S Pl'- 

{ai.ai^ak^ak^)^^ + 2wkMui)aiu2)\ < ©(t*^/^"^") (10.20) 

and for ui = 0,U2 £ Pl, 

{ai^al^ak,ak,)^^+2wkMuiHu2) < 0{q'''^-^'^ N). (10.21) 
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Furthermore, the smoothness of y implies |l/u^_fc^ — Vfe^ I <q^^^. So summing 
up ui,U2 ^ Bl{ui) and ki, k2, we obtain 

Ki-fci(4i42"fci«fc2)*c +2 a{ui)a{u2)\\Vw\\i 

ui,U2,k\,k2 ui,U2GPlUPo 

Y l^fcl^fcil^""^ - ll^^lll 



< 



_ 2 Y^ a{ui)a{u2 

U1,U2£PlUPo 

+ const. Y2 a{ui)a{u2) 

ni,n2:«2£-BL(wi) 



(10.22) 



One can see the first Hue of the right side is less than e^N'^ /2. Here is 
independent of a and lim^es = 0. With the bound a{u) < nic for u G Pl, 
we can obtain that the second line of the right side is also less than e^N"^ /2. 
So we arrive at the desired result (jlO.Sp . 

At last, we prove (jl0.4p concerning ui^2 ^ Pl, ui "^2 and /ci_2 G Ph- 
By definition of Ma and fa, if (/^lai^o^ au2ak2\l) ¥^ and f3,j Ma, then 
7 ^ Ma{ui) U Ma{k2), (3 i Ma{u2) U Ma{ki) and 



I /a (7) I < const. 
This implies 

/a(/?)/a(7)(/3|ati4i"«2afe2l7) < const. a(ni) 



1/0 



Wk2 



(10.23) 



|/„(/3)|2 (10.24) 



Summing up /3 ^ Ma(u2) U M„(/ci), with the bound on |/a(/?)P dSS]), 
we have 



('^Ii4i"«2afc2)*a < const. a(ni) 



^'fc2 



/3^Mc(«2)UM„(fci) 

< a(ni)|Vw^l^'^"^'' (10.25) 
At last, using \wp\ < 47ra|p|~^ and ~ |/c2|, we have 



E 



Ul,U2&Pj^,kl,k2&PH 



(4i4i^"2 0A:2)'I'o 



< ^a(ui)^>=^-9^A 

■ui 

< A2p5/2 



(10.26) 
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11 Proof of Lemma [TO 

In this section, we will prove Lemma [TO] involving interaction energy between 
particles with momenta in Ph- We will show that the only contribution to 
the accuracy we need comes from four high momentum particles, to be 
computed in Lemma [T6l (|11.3p . We start with separating {Hhh)^^ i^to the 
main terms and the error terms. 

Define Ma{ki, k2, k^, k4,ui,U2) C Ma Ma as the set of (/?, 7)'s where 
P and 7 can be created from the same a G Ma as follows, 

Ma{ki,k2,k3,k4„Ui,U2) (11.1) 



= {(/?, -f) e Ma0Ma:3ae Ma s.t. All'l'^a = p and All^a = 7}, 



where ki,k2,k^,ki £ Ph and ui,U2 G Pj^. For /ci , A;2 , ^3 , ^4 6 Ph, ui,U2 G 
P^, we define Au^^uiMMMM 



With (lll.2p . we can separate the expectation value of Hhh hito two parts, 
main term (Lemma I16p and error term (Lemma I17p . 

Lemma 16. Summing up ki, /c2, k^, k^ G Ph, h / kj for i / j, ui,U2 G Pj^, 
we have 



A. 





(11.2) 



MM,ksM-Na\M-'\\Vw''\\i <y^'A, (11-3) 



where £3 is independent of a and lim^^o £3 = 0. 



Lemma 17. Define Ma{ki,k2,k^,k/i) as follows, 



Ma{ki,k2,k3,k4, 



) = ^ui,u2ePiMa{ki, k2, ks, k4,ui,U2). 



(11.4) 



Then we have 



ki&Pn il3,-y)<^Ma(ki,k2,k3,k4) 





Here ki ^ kj for i ^ j. 
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11.1 Proof of Lemma 1101 

Proof. By the definition of Ma, we have f3{k) < 1 when k G Ph and f3 G M^. 
So the expectation value of a\_^a\^akgak4 must be zero when ki = k2 or 
ks = k^. Together with the definition of Hhh-, we can rewrite {Hhh)^^ ^.s 

{Hhh)^^= Y1 E Vk,^ks^-'U(3)7JT){p\alalak,ak,h) (11.6) 

On the other hand, if /3,7 G and {P\a\^a\,^ak'j,o,kA\l) / for some 
^1,2,3,4 G Pff, then by the fact PL{f3,a) = Pl{'^,oi) is non-trivial (Def. 
[5]), there at most one pair oi {ui,U2} satisfying 

(/?,7) G Maiki,k2,k3,k4„ui,U2) (11.7) 

So combine (jll.Sp and (jll.Sp . with jVfc^-fcgl < Vq, we obtain the desired 
result (f5l^ . 

■ 

11.2 Proof of Lemma [TBI 

Proof. We start with giving the bounds on Aui,u2,A;i,fe2,fc3,A:4- 

Lemma 18. When ui,U2 G Pl a'^c^ = ±^2 or U2 G Bl{ui), for any 
ki G Pff, we have 

■'^ui,U2,ki,k2,k3,k4 (-^-'-•^) 

In other cases, Ay^^^^^^^iMMM is hounded by (Pq = {0}) 

\Aui,U2M,k2M,k4 - a{ui)a{u2)Fa{ui,U2fwk^Wk.^A''^\ (11.9) 

{g^/^A~^a{ui)a{u2), ui,U2 G Pl 
g^/^A-^Na(u2), ui G Po,U2 G Pl 
gi/SA-^N^, ui = U2€Po, 

where Fa{ui,U2) = 1 when ui = U2 = 0, otherwise Fa{ui,U2) = 2. 

Proof. Proof of Lemma [TH] 

First we prove (jll.Sp . One can see that it follows the definition of 
A.ui,u2,ki,k2,k3,k4 and the result 2 in Lemma[Sl 

Then we prove (jll.Op when ui,U2 G Pl- If (|11.7p holds, by the def- 
inition of Ma {ki,k2,k^,k4,ui,U2) in (jll.ip . we can see that there exists 
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a £ Ma, -^kl'k2'^ ~ 1^ ^"^^ •^kl'kl'^ ~ 7- With definition of fa, when 
a E <^i=iMa{—ki), we have 



= -Fa{ui,U2)y/a{ui)a{u2)A ^ -Wk^^J -Wk^faia) (11.10) 

fa{l) = -Fa{ui,U2)y/a{ui)a{u2)A''^ ^-Wk^^-Wk^faia). 

And when a ^ <^i=iMa{—ki), we have the following bound on |/a(/3)/a(7)|, 

4 

\fa{P)fa{l)\ < MuiMu2)A-^ll\V^Mfa{a)f (H-H) 

1=1 



On the other hand, if ki G Ph for 1 < i < 4 and 

(3,-/ e Ma and (/3|4^ 4^ 0^30^417) / 0, 



(11.12) 



then by the definition of Ma, we have (3{ki) = P{k2) = 1 and 7(A;3) 
7(A:4) = 1. This impUes 



^Ii42afc3«fc4l7> = 1 



(11.13) 



Combine (Ill.lOp . (|ll.llll and pi.l3p . we obtain that if AU'l^a = P and 
•^kl'kA^ = 7) when a £ Df^^Mai-ki), then 



/Q(/3)/a(7)(/3|4i42"fc3«A:4l7> = -^a(^ii,'U2) a{ui)a{u2)A Y{^/^w^,\fa{a)\ 

i=l 

(11.14) 

When a ^ n^^;^MQ,(— fcj), using (|6.5p . we have 

XI /a(/3)/a(7)(/3|4i42"fc30fc4l7) < const. ^^/^a(ui)a(u2)A"^ 

(11.15) 

(11.16) 



So with (jll.l4p and (I11.15p . we can see 



A. 



ui,U2,ki,k2 



MM + 0{Q^/'^)a{ui)a{u2)A~ 



Fa{ui,U2fA Q:(ui)Q;(n2)|/( 



i=l 



Where ^ is defined as the set 

A = {aeMa: Al[%la = (3 ^ Ma, A''^^^'^ = e Ma, a e nUMa{-k.i)} 
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With the definition of Mq, denoting Ui, kj^s as pi,p2, ■ ■ ■ Pe, we obtain 

A = ri=iM^{Pi) n r\UiMa{-ki) n M^i'"^ (11.17) 

Here ^ '"^ is defined in Lemma [13] (j6.25p . Using the results in Lemma [T3l 
we have Yla&A \fi^)\'^ bounded by 

1< j;i/(a)P<l-0(f?'/') (11.18) 

On the other hand, using (jl0.13p . with the fact |A;i+A:2| = 1^3+^41 < q^^^q~^, 
one can bound the Yl^^i \/—Wk^ in (I11.16P as follows 



1=1 



< 0{g^/^-'^) (11.19) 



Inserting pi.l9|) and (I11.18D into pi.l6p . with the fact a G M„(ni), i.e 
a{ui) = a{ui){i = 1,2), we arrive at the desired result (jll.Op . 

Similarly, using the bounds on Qa{0) and Qa(0,0) in (j6.29p and (j6.30p . 
one can prove (jll.9p when one of Ui belongs to Pq or both of them belong 
to Pq. m 

With (jll.9p , summing up ki,ks,ui,U2, one can easily obtain the desired 
result ([TO]) . ■ 

11.3 Proof of Lemma [TTl 

Proof. As in ^16j, to estimate the error term of the interaction of particles 
with high momenta, we need to use a new tool. We start with defining the 
set Ma{a, s, {vi, ■ ■ ■ , vt}). Let vi, ■ ■ ■ ,vt £ Pl and being in different small 
boxes Bl, i.e., 

BL{v^) ^ BUvj), for i/j. (11.20) 
For non-negative integers s, t satisfying s + 1 £ 2N and a £ Ma, define 

{{s+t)/2 m 
i=m+l i=l 

(11.21) 

where the tij's and kiS satisfy 

1. Ui G Pr, Pi G Pff, 1 < i < s + t. And Uj = for i < 2m. 
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2. {ui, 1 < i < s + t} is a, permutation of s zeros and {vi, ■ ■ ■ , vt} 



3. for any fixed 2m + l < j < s + t, a G Ma{—pj), i.e., a{—pj) = a{—pj). 
And pj 7^ — Pi for 1 < i < s + t. 

. We note: if /?, 5 G Mq, have the form in (jll.2ip . then it is easy to check 
n « ^ A C {1, • • • (. + t)/2} (11.22) 

By this definition, if (jll.121) holds, then P{u) = j{u) for any u £ Pj^, so 
there at least exists one Ma{a, s, {vi, 1 < i <t}) satisfying 



P and 7 G Maia,s,{vi, - ■ ■ ,Vt}) 



(11.23) 



E.g. Using Lemma [H we can see (|11.23p holds when we choose a = a, 
{vi,--- ,vt} = PL{P,a) = PL{'y,a). 

Then, for any Ma{a,s,{vi, - ■ ■ ,ft}), we define N{a,s,{vi, - ■ ■ ,ft}) as 
the set of the pairs 7) such that 

1. /3, ^ £ Maia,s,{vi,- ■ ■ ,vt}) 

2. there exist ki,l < i < 4 satisfying (jll.l2p but 

(^,7) ^ Ma{ki,k2,h,h). 
Here Ma{ki, k2,k^,ki) is defined in (jll.4p 



3. for any other a' ,s' ,{v'i, - ■ ■ if/?, 7 G Ma{a' , s' ,{v'i 

then 

s + t<s' + t' 



(11.24) 

,</}), 
(11.25) 



We assume ()11.23p and ()11.12p holds. Clearly, s + t = 2ort = implies that 
(/3,7) G Mq(/ci, A;2, /ca, /C4). Hence if A/'(a, s, •• • ,ft}) is not an empty set 
then 

s + t > 4, and t > 1 (11.26) 

By definition of iV(a, s, {vi, • • • , ft}) and ()11.13p . we can bound the left 
side of (jll.Sp as follows {ki ^ kj for i ^ j) 



E ^oA-i /«(/?)/«(7)(/3|4^4^afc3afcj7) (11-27) 
VbA"-^ |iV(5, ,i;t})| max 

~ r 1 Pn(^Ma{a,s,{vi,- ,Vt}) 



< 



faW)fah) 
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where \N{a, s, {vi, ■ ■ ■ ,vt})\ is the number of elements in this set. By 
definition of fa and a{u) < nic = Q~'^^ for any u G Pl, we have: if 
/3,7 e Ma{a,s,{vi, - ■ ■ ,vt]) then 

|/a(/3)^(7)l <const.* 

Then with the fact \wp\ < 47ra|p|^^ and a(0) < A^, we obtain 

< const. *+^(^>i-2^)^(^>~5'')*|A|-Va(a)P (11-28) 
So, the right hand side of (|11.27p is bounded by 

a,s,{vi--vt} 

(11.29) 

Define N{a, s, t) and N{s, t) by 

A^(a,s,t)= max {\N{a, s, {vi, ■ ■ ■ ,vt})\} (11.30) 

{hi,--- jft} 

7V(s,t) = max{7V(5,s,t)} (11.31) 

a 

With A^(5,s,i) and N{s,t), we can bound (|11.29|) by 

dirSZD < (HOSI) <|;i/(5)P 5] iV(5,5,i)£.^(^^-6'')*+^|A|-*-i 

a,s,t {vi---vt} 

<E E A^(s,t)i?'(i?~''')*+1Ar*"' (11-32) 

■5,* {Vl--Vt} 

For fixed t, the total number of sets {vi ■ ■ -vt^Vi G Pl} is bounded by 
l<(A£,7?^3)*(t!)-i<(^i-3'')*|A|*(t!)-i 

{fi---ft} 

On the other hand t is bounded above by the number oi B^s in Pl, i.e., 
t < |PL|/max{|Bi|} < const, g^'^''-^''^ , 

i 

where \Pl\ and \B^j^\ are the volumes of Pl and the small box -B^'s. Combine 
with ()11.26p . we have 

(fTL27|) < ^ Yl N{s,t){Q^-^'^y+'\A\'\tl)'^ (11.33) 

t=l s:s+i>4 

We claim that N(s, t) is bounded with the following lemma, which will be 
proved in next subsection. 



a(0) 



max{|u;fc|}^+Va(«)l' 
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Lemma 19. For any N{a, s, {vi, ■ ■ ■ , Vt}), s + t > 4 and t > 1, we have 

\N{a,s,{vi,--- ,vt})\<tlt^T'>\A\'i'+\g-^Y+^ (11.34) 

Together with (jll.33p . we obtain 

(fTr27|) < E E (£''"'°'')'''*i^^^|A|'*' (11-35) 

t=l s:s+t>4 

^l-47,-3xr^ 

i=l s:s+t>4 

Then < 1/2 imphes q^''^^''A^/'^ < 1 and £.i-i°^t3/4^i/4 ^ therefore, 
we arrive at the desired result 

(fTL27D < 0(1) < ^^A (11.36) 

■ 

11.4 Proof of Lemma [T9l 

We now prove Lemma [T9l 

Proof. Since (/?, 7) S A^(a, s, {vi, • • • , ft}), we can express them as (|11.2ip : 

(s+t)/2 (s+t)/2 

P= U -4^-,^' 5, 7 = n S (11.37) 

1=1 i=l 

We note that for any 1 < i < {s + t)/2 

{q2i~i,q2i} ^ {ki,k2} and {q2i~i,q2i} 7^ {ksjh}, (11.38) 

otherwise one can see (/?, 7) E Ma{ki, k2, k^, k^), which contradicts with 
(/?,7) £ N {a,s,{ vi,--- ,vt})- 

From (|11.12p . the sets {^i,--- , q2s+2t} and {qi,--- , g2s+2t} are very 
close, i.e., 

{qi, ■■■ , q2s+2t} - {kl, k2} = {qi,--- , Q^2s+2t} - {A;3, k4} (11.39) 

Denote the common elements in {qi} and {qi} by pi, p2, Ps+t-2- 
Then we have 

{qi} = kl, k2, pi, P2, Ps+t-2, (11.40) 
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{qi} = k3, h, Pi, P2, Ps+t-2, (11-41) 

We now construct a graph with vertices {ki, k2, k^, k/i,pi,l < i < s + t — 2}. 
The edges of the graphs consisting of (3 edges {q2i-i,q2i)-, l<i<(s + t)/2 
and 7 edges (g2j_i, ?2j), 1 < ^ < (s + 0/2- From (jll.l2p . we know each 
ki{l < i < 'i) touches one edge and each pi{l < i < s + t — 2) touches two 
edges. So the graph can be decomposed into two chains and loops. Thus 
there exist /, mj G Z and < mi < m2 < < m/ = s + 1 such that 

chains [k^^P^^P^^P^- --P^rn.-^ ^ f2( or k,) 

y kz < > P2mi < > P2mi+1 ' ' ' P2m2-2 < > k^^ Or ^2) 



loops < 



P2m2-1 ^ > P2m2^ ^P2m2+1 ' ' ' P2(m3)-2 < ^ P2m2-l 



L P2mi_j-1 < > P2m,_i < > P2»n,_i+1 " ' ' P2(mi)-2 ' ' P2m,_i-1 

Here we have relabeled the indices of p and do not distinguish /? edges and 7 
edges. We also disregard the obvious symmetry ki k2 and k^ ^ k^. Due 
to the condition (jll.25p and the facts a) = Pl^I, a) is non-trivial(Def. 

[5]), the length of the loop must be 4 or more, i.e., each loop has at least 4 
edges and 4 vertices, i.e, 

rui-i + 2<mi for 3<i<l (11.43) 

The inequality (jll.38p implies 771-2 ^ 2. Together with mi = (s + t)/2 and 
pi.43p . we obtain 

I < {s + t)/A + l, t>l. (11.44) 

Without loss of generality, we assume mi — mj_i is creasing with i > 3, i.e., 

for 3 < i < j < I 

mi — mi-i < mj — mj-i (11.45) 

Denote by N{a,s,{vi, - ■ ■ ,vt},l,{mi, - ■ ■ ,mi}) the set of all pairs (/?, 7) 
having the graph above and we now estimate the number of elements of this 
set. 

Using the notions Wi = {w2i-i,W2i) and Wi = {'W2i~i,W2i), we can add 
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the information between fcj's and pj's into the graph as fohows 

ki < — > pi i — > p2 < — > P3 ■ • •P2mi-i — ' hi or k2) (11.46) 

k3< > P2mi < > P2mi+1 ■ ■ ■ P2m2-2 ^ > k2[ Or ki) 

Wm,2 + 1 Wm,2 + 1 W^mg 

■ ' "2m2 ^ > P2m2+1 ' ' ' P2(ms)-2 ' ' P2m2-l 



P2m,_i-1 < > P2m,_i < > P2m,_i+1 • • • P2(m,)-2 ^ ' P2mi_^-1 , 

where Wi^s are the union of s zero's and {vi,--- ,vt}, so are w^s. More 

specifically, if A •s-^ B appears in the graph and W = {C,D), then the 
operator appears in (|11.37|) . So we have 

A 3-^ B <^ A + B = W2i-\ + W2i, 

so as VF's. With this relation, we can see that /? and 7 is uniquely determined 
by the structure of the graph, Wj's, Wj's and one ki or pi for each loop or 
chain. 

To bound |iV(3, s, {ui, • • • , ft}, /, {mi, • • • , m^})!, we note that the sum 
of momentum (pi's) in each loop is zero. Thus we can count the number of 
graphs as follows. 

1. choose the positions of zeros in /3 edges. The total number of choices 
is less than 2*+**. 

2. choose the positions oiv\ - ■ ■ Vt'vcL (i edges. The total number of choices 
is t\. 

3. choose the positions of zeros in 7 edges. The total number of choices 
is less than 2*^^^ again. 

4. choose the positions oiv\---vt in 7 edges. We call a loop trivial if all 
the momenta associated with 7 edges are zero. The number of trivial 
loops is at most s/4 since there are at least two 7 edges per loop. 
Hence the number of non-trivial loops is at least / — s/4. Thus we only 
have to fix v in at most t — {I — s/i) edges and the number of choices 
is at most t^-^+^/^. 
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Thus we obtain 

\N{a,s,{vi, - ■ ■ ,vt},l,{mi, - ■ ■ ,m/})| (11.47) 

< (const. )*+^t!t(*+'^/4-') {g-'^'^Ay 

< (const. )*+^t!t(3*/4) (^"3„^)*/4+./4+l 

where we have used (|11.44p Since 

\N{a,s,{vi,--- ,vt})\ = ^ \N{a,s,{vi,--- ,vt},l,{mi,--- ,mi})\ 

I {mi,---,mi} 

and 

^ ^ 1 < const. (11.48) 

/ {mi,---,m;} 

we have proved ()11.34p . ■ 

12 Proofs of Lemmas 1, 2, 3 
12.1 Proof of Lemma [T] 

The proof of Lemma [1] is standard and only a sketch will be given. We first 
construct an isometry between functions with periodic boundary condition 
in A = [0,L]'^ and functions with Dirichlet boundary condition in A* = 
[-£, L+e]^, where L = g'^^/^'^ and £ = g''^^/^'^^. We note, by the definition 
of g* in (I33D . 

\Mq = \^*\q* (12.1) 

Denote the coordinates of x by x ^ '). Let /i(x) supported 

on [—£, L + £]^ be the function /i(x) = q(x^^^)q{x^'^^)q{x^^'^) where 



(x) 



' cos[(x-^)^/4£], \x\ <£ 

1, £<x<L-e , 

cos[ix-iL-£))TT/4£], \x-L\<£ ^ ' 

0, otherwise 



The function q{x) is symmetric w.r.t x = L/2. Due to the property of cosine, 
for any function (j) with the period L we have 

/ |/i0(x)|2dx= / |</)(x)|^dx (12.3) 
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Thus the map </> — > h(f) is an isometry: 

-^Periodic ([0) ~^ -^Dirichlet ([~^' ^ + ^]^) ■ 

Let x(x) be the characteristic function of the ^-boundary of [—£, L + £]^, i.e., 
x(x) = 1 if < i for some a = 1,2 or 3 where is the distance 

on the torus. Then standard methods yield the following estimate on the 
kinetic energy of hep 

[ |V(M(x)|2 (12.4) 

Jxe[-£,L+i]^ 

< [ |V,/.(x)|2 + const, /x(x)|0(x)|2 

ixe[o,Lp J 

The generahzation of this isometry to higher dimensions is straightfor- 
ward. Suppose ^'(xi,--- ,xjv) e Aperiodic ([0'^]^^)- Here 

N = \A\g = \A*\g* (12.5) 

Then for any u S M^, the map 

N 

r'i^) :=^'(xi,--- ,Xi,--- ,XN)Ylh{xi + u)h{yi+u)h{zi+u) (12.6) 

i=l 

is an isometry from Ll^^,^^,^ ([0, Lf^) to Ll,^,^^,^^ (^[-l-n,L + £- n]^^) . 
Here the x's in ^ are coordinates on the torus. Clearly, has the property 
([123]). 

The potential V can be extended to be periodic by defining V^{x — y) = 
V{[x — y\p) where [x — y\p is the difference of x and y as elements on the 
torus [0, L\. Since V is nonnegative and has fast decay in the position space, 
we have V{x — y)< V^{x — y). From the definition of J-^, we conclude that 

» N N 

/ ir'(^)iv(xi-x2)nrfx, < / iM'|2F^(xi-x2)n^x, 

J[-t.-u,L+l-u]^'^ fj[ J[0,L]3iv fj-^ 

Therefore, the energy of two boundary conditions are related by 

N 

< (Fiv)^ + const. r2^(x(xi + n))^ (12.7) 
1=1 
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We note is operator on pure state ^ . It can be generalized to operator 
Q'^ on state V as follows. For any state of N particles in [0,L]^ with 
periodic boundary condition, we define 

g"(r^) :=^r^(^)t (12.8) 

So F^ = ^"(F^) is a state of N particles in [— ^ — n, L+£ — n]^ with Dirichlet 
boundary condition. With ()12.ip . one can see 

G^: F^(^, A, /3)^F^(^*, A*, /3) (12.9) 

Using (|12.7p . we have: 

N 

TV g"(F^) < Trfl-TvF^ + const, r ^ ^ Trx(xi + n)F^ (12.10) 

i=\ 

Averaging over u G [0, L]^, we have 

L"^y" (TrF7vg"(F^)) (in <TrF7vF^ + const. r^L^^iV (12.11) 

So for any F^ there exists at least one u such that 

TrFiva"(F^) < TrFivF^ + const. ^{^) (12.12) 

On the other hand, the fact .7-""( (|12.6p ) is a isometry implies that ^"(F^) 
and F^ have the same von-Neumann entropy, i.e., 

5(g"(F^)) = 5(F^) (12.13) 

Combine (jl2.12p and (112. 13p . we obtain A/ the free energy difference 
between ^"(F^) and F^ is less than const. A/^(;^). With the choice L = 
^-41/60 g^j^^ ^ _ ^-41/120^ ^j^g error term is negligible to the accuracy we 
need in proving Lemma [H This concludes the proof of Lemma [H 

12.2 Proof of Lemma El 

It is not easy to construct Fq (the state of A'' particles) directly. So we 
start with constructing a state Fjf in Fock space. Then pick up the useful 
component of Fjf and revise it to Fq . 
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First, let Bjr be the standard basis of the Fock space ^(A) as foUows 



a):\a) = Co, J] (4r(^>), a(A;) G N U {0} [> , (12.14) 



where Cq, is a positive normahzation constant. We define a new Hamiltonian 
H of free particles in [0,L]'^ with the following revised Bose statistics, i.e., 

1. The number of the particles in single particle state \k) is nonzero only 
when k £ PjU Pl- 

2. The number of the particles in single particle state |A;), k G PlL) Pj, 
must be no more than Ck, which will be chosen later. 

With the definition of fi in (j2.8p . we define Fjf as the grand-canonical 
Gibbs state of H with the chemical potential fi{g, /3) < and temperature 
where 

^^^(l_L-i/2) = ^(l_o(^i/3)) (12.15) 
and Cfc is chosen as follows (Recall rric = Q~^^) 

Ck = { ^ ^ , (12.16) 

y nic k £ Pl 

where E^^^ is defined as k'^ — fi{Q,(3). We note that /? ~ g~'^^^ implies 
PEk,^Ck > |log^|. 

We claim that the state Fjp has the following properties: 

Lemma 20. The free energy per volume ofTjr is bounded above by 

fir:F)<foiQ,m-o{g'/')) (12.17) 
In most cases, the total particle number ofVjr is less than N = gA, i.e., 

N 



Tr^^F^ > 1 - ^ (12.18) 

m=l 

Here F^ are the components ofVjr on TLm, i-^-, 

oo 

r^ = J2 : ^ Hm (12.19) 



m=0 
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Similarly, in most cases, the total particle number of Tjr is very close to 
mm{g, gc}A. We have 

E Tr^^r^ >^-e (12.20) 

\m-min{g,gc}A\<N g''-/'-^ 

Proof, proof of Leml20] 

First, we prove ()12.17p . by the definition, the free energy of Tjr is 



-1 



2^ lo! 



(12.21) 



k^PLUPi \ kePLUPj 

With the definition of Pj and Pl, one can remove the Ck terms, add the 
k ^ Pi yj Pl terms and check that (|12.2ip is equal to 

\ fce{^)3,fe7^o ^ ^ fce(2|Z)3,fc^o j 

(12.22) 

Then with the choice L = q"^"^/^^ and the definition of free energy /o in 
([2:6]) and (HZD, we have 

([12:221) = Uq, /5)A(1 + oio"^^)) (12.23) 

Combine with g = q{\ + o{q^/^)), we obtain the desired result (|12.17p . 

Then we prove ()12.18p . Let n[k) denote the number of the particles in 
one-particle-state \k). So n{k) the average of n{k) is equal to Tra|.a^, Fjf. 
By the definition, the average total number of particles of Vjr is equal to 

y ^= y —1 y \+^^^ — (12.24) 

k&PlUPL k&PiUPL kePLUPi 

Similarly, with L = ^-^^eo ^nd (iEu^^Ck » | log q\, one can easily prove: 

(112:21) = min{^,^JA(l + 0(^>-i/3L-Mog^)) (12.25) 

On the other hand, since n{k) are independent random variables and they 
are bounded in (|12.16p . we can use Hoeffding's inequality [6J to estimate 
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the distribution of the total particle number of Tjr. With n{k) < Ck and 
Hoeffding's inequality [6], we obtain that the probability of find more than 
particles in Tjr is bounded above by 



P ^ n{k) > n \ < exp < 



By the definition of (|12.16|) . we have: 

k&PiVjPL 



(12.26) 



(12.27) 



On the other hand, with the fact g — q = qL and ()12.25p . we can see 
that 



(12.28) 



Inserting L = q~'^'^/'^° , (|12.27p and (|12.28p into (|12.26p . we obtain the desired 
result (|12.18p . And (112. 2Up can proved similarly with p2.25p and p2.27p . 



By Lemma [20l there exists rriQ s.t. 

mo < iV and \mo - m.m{g, Qc}A\ < g^^^N (12.29) 

and the free energy of T^" is less than fo{g, P)A{1 — o{g^/^)). 

Then adding — ttiq {N = gA) particles with momentum zero into the 
system described by F^", we obtain a new state Tq of free particles. The 
state Tq always has — mg particles with momentum zero and the free 
energy of Tq is also less than /o(£>, /3)A(1 — o{g^^^)), i.e., 



A-i < o{g') 

Furthermore, by the definition of Tj^, Tq has the form: 



Tr-AFo + -5(Fo)-/o (£»,/?) 



^0 = 'Y gaig,P)\a){a\, a{0) = N - niQ and = 1 



(12.30) 



(12.31) 



We note: if a{k) > Ck for some k £ Pj U Pl, then gaio^P) = 0. This 
property implies that the total number of the particles with momentum in 
Pj is o{N), i.e., for any a £ M 



aeM kePi 



(12.32) 
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Together with (|12.29p and the fact a{k) < rric for a G Pl, we obtain (j3.23p . 
At last we prove ()3.22p . First with the structure of Fq, we have 



TrW]v,A^ro= J2 9a{g,P){a\V\a) 



(12.33) 



= ^ga{Q,f3)i Yl VoK-\a{kf - a{k)) 

a&M ^ feePoUP/UPL 

+ 2{yo + Vk-k')^-^a{k)a{k')\ 

k,k'ePoUPiUPL ^ 

Using the smoothness of V and <C 1, we can replace Vh-w with 

Vo without changing the leading term. Then with the cutoff C^'s, the fact 
Q!(0) = iV - m and (|12.29p . we have 



lim \lWT^\Q-'k~^ = Vo{2 - [1 - R[f5] ]^) 



(12.34) 



Combine with (jl2.30p . we obtain ([322]). 
12.3 Proof of Lemma [3] 

Proof. Since the states |a)'s in M are orthonormal, we can rewrite the en- 
tropy of Fq in lemma [2] as 



5(Fo) = - ^ ga^ogga 



(12.35) 



For S{T), we define A^o as Aoo = l^")(^"l|loo rewrite F as 

ya- 



With the fact TrF = 1, i.e., ^ (7^ = 1, we have 
5(F) = - log Aoo- Aoo Tr 



(12.36) 



(12.37) 

With Berezin-Lieb inequality [11], [2], we obtain 

5(F) > - log - <?a log ga = - log Aoo + 5(Fo) (12.38) 

On the other hand, we claim the following lemma 
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Lemma 21. 



lim loe 



(12.39) 



Insert this lemma into (|12.38p , we arrive at the desired result (j3.27p . 
12.3.1 Proof of Lemma [2T] 

Proof. With the fact: for any hermitian matrix M = Mij, 



'■i3\ ( ^ 



\\M\\oo < max |E 
we can bound || YlaeM l^a)(^«|||oo as follows (Recall M in Def 2.) 



E l^a)(^«|||oo < max< 



/3eM 



E El^/5|*")^*-I^)l[ (12.40) 



7GM 



< maxJ E K/?l^«)l[-maxi E 



With the fact ^q, is the linear combination of states in Ma C and 
l^'rv) are normalized, we claim 



log (max I E ) < g'-'^-^^^ 



(12.41) 



log I max \ E K7l^a)l ^ ) < g''^^-^^^ + q-^^-^^^ (12.42) 



,7eM 



First, wejH-ove (I12.4ip . We know K/Jj'I'a)! t^J) implies K/3|^'a)| ^1, a G M 
and P £ Ma- With the definition of M and M^, if a € M, /? G M^, then 



/3(n) = q(u) for u £ Pj 
P{u) < a{u) for u e Pl 
a{u) = for ti G 



(12.43) 
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and for any fixed small box Bl{i = 1,2,...) in Pl, P{u) is very close to 
a{u), s.t. 

^ \Piu)-aiu)\<l (12.44) 

Now let's count, for fixed /?, how many a € M satisfying /? S Mq. This 
number must be less than the a's satisfying (jl2.43p and (jl2.44p . By the def- 
inition of -Bl's, the total number of B^s is less than const, g^^'^'^^^^^ . And 
for any B\^, \B\^\ the number of the elements in B\^ is less than const. g'^^^A. 
Therefore, for fix /? E M, the total number of a G M satisfying (3 G Mq, is 
less than 

(const.^.^^-^A)™"^*-^ " (12.45) 

Together with the fact < 1, we proved (I12.41|) . 

Then we prove (jl2.42p . Similarly, using the rule 2 of Def. 3, we can 
count, for fix a G M, the total number of 7 G M, s.t. |(7|^a)| / is less 
than 

(const. £,3-. A)--*- (const. a)™-*- ^"^^-^^^ , (12.46) 

which implies (|12.42p . Inserting (|12.4ip and (|12.42p into (|12.40p . we obtain 
the desired result (|12.39p . ■ 
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